Population Dynamics in Cold Gases Resulting from the Long-Range Dipole-Dipole 

Interaction 
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We consider the effect of the long range dipole-dipole interaction on the excitation exchange 
dynamics of cold two-level atomic gase in the conditions where the size of the atomic cloud is large 
as compared to the wavelength of the dipole transition. We show that this interaction results in 
population redistribution across the atomic cloud and in specific spectra of the spontaneous photons 
emitted at different angles with respect to the direction of atomic polarization. 
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I. INTRODUCTION 

Cold atomic gases give a possibility to address collec- 
tive atomic quantum states, where an elementary atomic 
excitation is coherently distributed among a large num- 
ber of individual atoms. Such systems have been consid- 
ered in the context of Quantum Information coher- 
ence protection [5], individual photon manipulations fsj 
etc. For cold Rydberg gases that have large dipole mo- 
ments, one has to allow for the effects associated with 
strong dipole-dipole interaction among atoms 
Considering the regime typical of current experimental 
setting, one usually takes into account the regular static 
dipole-dipole interaction V ^ among pairs of the 

atoms, which from the view point of Quantum Electro- 
dynamics is a result of exchange by a virtual strongly 
off-resonant vacuum photon with a typical wavelength of 
the order of interatomic distance R . This sort of in- 
teratomic coupling can be attributed neither to the short 
range interaction, nor to the long range ones, since the 
integral V{R)R^dR corresponding to the average bi- 
nary interaction logarithmically diverges at both upper 
and lower limits. This circumstance results in a number 
of interesting dynamic phenomena [1| , such as incomplete 
decay of a single atom population at the limit of long 
times, which resembles the Anderson localization effect. 
In the same time, for the description one cannot make 
use of the continuous media model and the mean field 
approximation, which are applicable only for the long 
range interactions. 

Among experiments with Rydberg atoms, there are 
some [H that have been done in the regime where the 
size of a cloud of cold Rydberg atoms is of the order or 
larger than the wavelength of the resonant dipole active 
transition. In such conditions, the role of the radiation 
trapping, super-radiance [9|, and the dynamic long-range 
dipole-dipole interactions (~ 1/^) is dominant, and the 
continuous media model becomes applicable. In other 
words quantum states of the atomic ensemble and that 
of the resonant radiation field get strongly entangled, 
whereas the eigenstates of the compound atoms-|-field 
system correspond to atomic states dressed in the reso- 



nant radiation field. Driven by the curiosity we will now 
consider the atomic population dynamics in the simplest 
version of such a process. We make use of the mean-field 
approximation, and concentrate on an exactly soluble 
case of an infinite uniform continuous two-level media, 
ignoring the contribution of the static dipole-dipole in- 
teraction not conforming the requirements of the mean 
field model. 

We consider a static continuous media corresponding 
to the gas of Af two-level particles with the transition 
frequency w and polarization d along z-direction that 
are uniformly scattered in space at fixed positions. Since 
the static dipole-dipole interaction between the atoms 
can be viewed as a process intermediated by virtual pho- 
ton exchange, the effective Hamiltonian that governs the 
evolution of the combined system in the rotating wave 
approximation, is 
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are the collective operators of the exciton k, Uk = 
sin a is the coupling and a is the angle between 
the direction of the polarization d and the normal to 
the photon polarization plane. The coupling strength 
U\/k = dV ATihckn is a product of the photon vacuum 
field strength E ~ Anhck /V in a volume V —N/n 
and the collective atomic dipole moment d\07, where d 
is the atomic transition dipole moment, c is the speed of 
light, and k — |k|. 

We assume that at time t = 0, all but one particles 
are in the ground state. The consideration is equally 
applicable to the case where a group of particles locating 
in a volume with a typical size a, small as compared to the 
resonant transition wavelength 2ttc/u!, is in a coherent 
superposition of the individual excited states, such that 
the total number of excitations is one. In order to take 
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advantage of the uniform distribution, we take states 
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corresponding to a given wavevector k as the basis set in 
the coUective Hilbert space. Here |0) denotes the vacuum 
where ah particles are in the ground state, and the con- 
dition that the average distance rT^I'^ among the neigh- 
boring particles corresponding to the density n is much 
shorter than the transition wavelength l-ncjuj is implicit. 

Each state, ([3]), interacts only with the photon of the 
same wavevector, and the amplitudes '^■^ and (/Jk of the 
particle-|- field compound state |k)^ = i/^k |k, 0) -|-(^k |0, k) 
(where the second quantum number corresponds to pho- 
tons) satisfy the Schrodinger equation 

ifti/'k = C^\/fcsin(a) (^k H — 3^(i) (4) 

(27r) 



i?ii^k = ^{ck — tj)iy9k + U ^fk sin (a) tpi^ 



where Dirac delta fmiction 6{t) stands for the initial con- 
dition corresponding to the excitation location in the ori- 
gin, that is at the point r = 0. 
The exact solution of Q and ([5|) 
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with the Rabi frequency n = \J (^^^)^ + k (^)^ sin^ a 
suggests 

V'(r) = I /;2g-ifcrcose+it 



<p(r) = 
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for the coordinate dependent amplitudes, with d7 = 
dfcdr, where dF denotes the integration over solid an- 
gle and the integration over k starts at the point and 
goes to +(X). 

Long range dipole-dipole interaction corresponds to 
the exchange by photons close to the resonance domain 
/c ~ w/c. In this domain, the integrands ([H])-® have two 
branching points 
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in the complex plane of k where i7 assumes zero value. 
Therefore the contribution of the long range interaction is 
given by an integral along a contour presented in figure [1] 
circumventing these points, which can be expressed in 
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FIG. 1: Transformation of the integration contour for the 
inverse Fourier transformation integrals Q-Q. The initial 
contour Ci can be moved to the lower part of the complex 
plane of k. The part of the contour C2 circumventing the 
branching points (|10p accounts for the contribution of the 
long-range dipole-dipole interaction, while the remaining part 
of the contour accounts for the regular static dipole-dipole 



(5) interaction ~ 1/R 



terms of the Bessel functions J„ (x) . In order to shorten 
the notations we employ the units where w /c = \ hence- 
forth. In the limit L/ ^ 1 of a coupling, which is weak 
as compared to the energy of atomic transition quantum, 
the integration over dfc yields 

fiL We)eM«-.-e) ^ (11) 
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Ji (2U^^re(t^re)) e-'('---'^') , (12) 
7r^'yU^QH{re)QH{t-re)y. (13) 
J^{2U^^Te{t-re)) c-'<-^-'"") , (14) 



up to terms of higher orders in U. Here Qh (x) is the 
Heaviside-0 functions, while re and Ua stand for rcosG 
and U sin a, respectively. 

Integration over the solid angle in (|lip - (fT51) has to be 
done numerically in the polar coordinates associated ei- 
ther with the direction of the polarization d (angles a and 
(f)) or with the direction of the radius vector r (angles Q 
and $) as shown in figure [21 It has to be performed over 
the intervals {0, tt} and {0, 27r}, respectively. Angles a 
and 8 do not belong to a same reference frame and there- 
fore one of these angles has to be expressed in the frame, 
to which belongs the other one. For the purpose, one can 
make use of the relations 

sin^ a — 1 ~ (cos 9 cos 8 + sin 9 sin 8 cos $)^ (15) 



cos 8 — cos 9 cos a + sin 9 sin a cos ( 



(16) 



that follow from the invariance of the scalar product (kr) . 

Alternatively one may employ the cylindrical set of 
coordinates z — rsin^, p — rcos9, and 0, associated 
with the direction of the polarization d. 
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FIG. 3: Probability of the excitation w{r) = \4'{r)f for the 
infinite media and for times T = tU^ <^ 1 given by Ea. l|18p . 



FIG. 2: Two sets of sperical coordinates in the space of the 
vavevectors k: the coordinates {k, a, 0} associated with the 
direction of polarization d, and the coordinates {k, O, <E>} as- 
sociated with the radius vector r. For a given angle 9 betweeen 
the direction of the polarization and the radius vector, among 
these coordinates there exist relations (|15|) - (|16| ). The photon 
polarization E is on the plane perpendicular to k. 



II. POPULATION DYNAMICS AT THE 
WAVELENGTH SCALE 

We note that there are two size scales in the problem: 
the wavelength 27r (that is 27r/fc in the dimensional units) 
and a distance 1/U (that is c/U in the dimensional units) 
at which the light propagates during the Rabi period at 
resonance. Population distribution at these scales have 
to be analyzed separately. 

Let us first consider the wavelength scale. For short 
times t ^ the Bessel function can be cast in Taylor 

series and the first order yields 
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for the amplitude in . With the help of equation 
one can calculate the integral exactly and find 



V'(r) = 



tU^ (3r2 - 2e-''' (r^ + ir + l) 
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The distribution of figure [3] persists as long as tU"^ re- 
mains a small number. Its structure can be easily un- 
derstood when one notice that according to the second 
equation ([7]), the amplitudes of emitted photons are lo- 
cated in the domain of resonance where the the detunings 
(1 — fc) are small. These photons, after being emitted, re- 
sult in discarding of the corresponding spacial harmonics 
from the initial (5-like excitation profile '(/'(r) containing 
all harmonics with equal amplitudes. But such a dis- 
carding implies establishing of the population distribu- 
tion given by the discarded harmonics taken with the 
opposite signs. As the result, the sum of the (5-like dis- 
tribution and the latter distributions of ■0(r) no longer 
contains harmonics strongly coupled to the cooperative 
resonant radiation. In course the of time the harmonics 
with higher detuning are getting involved in the process, 
and the net population of figure [3] increases. For an ini- 
tial distribution of a small but finite size a, the relative 
amount of this population is given by the ratio of the 
length U'^t of the /c-interval of emitted photons and the 
length ~ 1/a of the fc- interval occupied by the initial 
distribution. 

In course of time t > U^^, the distribution |'0(r)|^ in 
the domain r 1 gets modified, as shown in figure 31 

One can find an analytic expression for the asymptotic 
form oi iP{y) at t — > oo. The main contribution to the 
integral comes from the domain of small cos0, that 
results from the exchange by the resonant photons with 
k almost orthogonal to r. Integration over dP can be 
carried out in the reference system associated with r with 
(18) the allowance for the approximate relation sin^ a ~ 1 — 



which results in the population |'0(r)| quadratically in- 
creasing in time, shown in figure [31 Note, that at long 
distances the probability corresponding to decreases 



as r~ , which means that the integral of the population 
over the volume logarithmically diverges at the upper 
limit being truncated just at the radius r = t. In other 
words, some part of the excitation is getting transferred 
by the resonant photons at a large distance r ~ t. 



(sin 9 cos $) and yields 
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One sees that the population distribution decreases with 
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photons detected at a distance much larger as compared 
to the typical size of the atomic cloud. 

The calculations relevant to the case of large distances 
rely on the Fourier transformed ^ and ^ 



sV'k = U Vk sin (a) (^k + 
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(k — l)ipk + uVksin (a) i/'k, 



while the resulting amplitude 
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(23) 



to have a photon with the wavevector k at an angle a to 
the direction of the polarization has to be averaged over 
an interval around fc = 1 of a width ^ 2t: / L correspond- 
ing to the uncertainty resulting from the finite size L of 
the system. The averaging has to be performed with a 
weight function n(K) given by the Fourier transformed 
particle density profile 



'n(K)U sin a dn 



87r3 (-£2 + + c/2 sin^ a) 
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For an analytic profile «(«;) this yields the spectral inten- 
sity 
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where the frequency deviation from the resonance 5lj ~ 
ejU sin a is scaled by the interaction. 

In figure [7| we present the probability of the photon 
detection as a function of frequency calculated for the 
particle density distribution n(r) = 1/ cosh^(|r| jV) with 
L = 4. This distribution stands as a physical approxima- 
tion to the delta function spatial distribution of excited 
atoms we have assumed throughout the derivation. One 
sees, that though the cooperative coupling between the 
particles and the electromagnetic field does not govern 
the population distribution over the ensemble of 

two-level particles of a finite size L < c/C/, it manifests 
itself in the spectrum of emitted photons that can be 
registered at distances large as compared to L. 



IV. CONCLUSIONS 



the ground state, has been excited in a small domain 
close to the center. In course of time the excitations gets 
redistributed over all the volume of the cloud due to the 
long range dipole-dipole interaction. The typical rate of 
the process corresponds to the l]^ — > AiKpn/h. For short 
times, the population distribution is given by (|18p and 
for tU^ ^ 1 it takes the asymptotic form (fTO|) . The pro- 
cess is associated with the cooperative spontaneous emis- 
sion of photons. The typical scale corresponding to the 
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FIG. 8: The long range dipole-dipole interaction results in the 
redistribution of the probability amplitude to be in the excited 
state, which occurs with the typical rate . The population 
distribution for small times and for long times are shown by 
the contour plots in the upper and the lower parts inside the 
circle representing the atomic cloud. In 3D the distribution 
has cylindrical symmetry around the polarization axis d. The 
process is associated with the cooperative emission of photons, 
and the typical scale of the spectrum of these photons is of 
the order of U . 



spectrum of emitted photons is given by the interaction 
strength U d^jAmmilh, and at a large distance from 
the atomic cloud it has a two-peak structure of equation 
(1^ . As a general conclusion, one may state that the 
presence of two essentially different energy scales {U and 
U^) can be considered as typical manifestation of collec- 
tive atomic phenomena. 



We conclude by presenting the overall picture of the 
process shown in figure [51 An atomic cloud, initially in 
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